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I discuss the von Neumann entanglement entropy in two-dimensional quantum Lifshitz criical point, 

namely in Rokhsar-Kivelson type critical wavef unctions. I follow the approach proposed by B. Hsu et al. 

[Phys. Rev. B 79, 115421 (2009)], but point out a subtle problem concerning compactification of replica 

boson fields: although one can define a set of new boson fields by linear combinations of the original fields, 

the new fields are not compactified independently. In order to systematically study boundary conformal 

field theory of multicomponent free bosons, I employ a geometric formulation based on compactification 

lattices. The result from the boundary conformal field theory agrees exactly with alternative calculations 

by J.-M. Stephan et al. [Phys. Rev. B 80, 184421 (2009)], confirming its universality as argued originally 

by B. Hsu et al. 
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£ . 1 Introduction 

o: 

Recently, stimulated by developments in quantum information theory, characterization of quantum many- 
body systems in terms of quantum entanglement is studied vigorously. Entanglement entropy is a quan- 
titative measure of the quantum entanglement. A natural problem then, is to divide the system into two 
^ ' regions and study the entanglement entropy between them. In general, short-range correlations, which are 

non-universal, give a contribution proportional to the area of the boundary between the two regions to the 
entanglement entropy. This is often called as the area-law contribution. In one dimension, the boundary 
consists just of points, and the area-law contribution is constant with respect to the system size or to the 
size of the regions. 

Calabrese and Cardy pQ studied the entanglement entropy between a region of length L and the rest of the 
system, in an infinitely long one-dimensional system. They demonstrated that, in a critical one-dimensional 
system described by a conformal field theory (CFT) with central charge c, the entanglement entropy has a 
universal contribution (c/3) logL. In addition to being interesting in itself, this has proved practically useful 
in determining the central charge of the system from numerical simulations. 

In two (and higher) dimensions, the non-universal area law contribution diverges and is the leading 
contribution to the entanglement entropy, as the size of the region increases. Nevertheless, there can be 
universal contributions as subleading terms in the entanglement entropy. Kitaev and Preskill [2] , and Levin 
and Wen [3] found that, in two-dimensional topologically ordered phases, there is a universal constant term 
reflecting the topological order, in addition to the non-universal area law contribution. By considering a set 
of geometries, one can cancel out the area law contribution to obtain the universal constant. This is also 
confirmed numerically. 

Classification of critical points in two-dimensional quantum systems is not as well understood as in one 
dimension. Nevertheless, there is a class of two-dimensional quantum critical points, for which the universal 
critical phenomena can be described precisely. 
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Some quantum systems in two dimensions can be related to two-dimensional classical statistical systems, 
by identifying the wavefunction in the quantum system with the statistical probability of the correspond- 
ing configuration in the classical system. This is a generalization of the Rokhsar-Kivelson wavefunction 
introduced for quantum dimer models. [I] 

When the classical system is at a critical point, the quantum system is also at a quantum critical point. 
A critical point of two-dimensional classical statistical systems is often described by a two-dimensional CFT. 
If this is the case, the corresponding quantum critical point of the quantum system is also governed by the 
same CFT. Let us call such quantum critical points as two-dimensional conformal critical points, following 
Ref. [5J. Although such a wavefunction is rather special, I can expect that a conformal critical point represents 
a certain universality class. An important subclass of such two-dimensional conformal criticalities is those 
correspond to the free boson CFT. Following Ref. [BJ, I call this as quantum Lifshitz universality class. 
Quantum Lifshitz universality class is in fact a one-parameter family of universality class, since the free 
boson field theory is characterized by a free parameter (compactification radius). 

Fradkin and Moore [5J, and subsequently Hsu, Mulligan, Fradkin, and Kim [BJ studied the entanglement 
entropy in two-dimensional conformal quantum critical points using replica trick and boundary CFT. In 
Ref. [5J, it was argued that the entanglement (von Neumann) entropy has the form 

S B = al-^(A X )\n(^j+0(l), (1) 

when the boundary T is smooth. Here c is the central charge of the CFT, Ax is the change in the Euler 
characteristics by the partition of the system, I is the length of the boundary T between the regions A and 
B, a is the non- universal coefficient of the area law contribution, and a is the ultraviolet cutoff. 

When r is smooth and Ax = 0, the logarithmic term vanishes. In Ref. [BJ, it was argued that in such a 
circumstance, the 0(1) term in eq. ([T]) contains a universal constant, similarly to the case of topologically 
ordered phases discussed in Refs. [21 [3] . 

Following these developments, Stephan, Furukawa, Misguich, and Pasquier [7] studied the same constant 
term in eq. (Q} using different analytical and numerical methods. They agreed that the constant is a universal 
quantity determined by the underlying CFT. In particular, the universality is confirmed for lattice models 
with different microscopic parameters. However, the universal constant obtained in Ref. [7j for the quantum 
Lifshitz universality class disagreed with the original prediction in Ref. [BJ . The disagreement is potentially 
serious, since if both derivations stand valid, it would imply a breakdown of the universality of the constant 
term in the entanglement entropy. 

In this paper, I aim to resolve the issue by re-examining the derivation of the entanglement entropy in 
Refs. [5j [BJ. I point out there are subtle problems in "changing the basis" technique employed to derive a 
fundamental formula in Refs. 0(6]. (I note that, although the fundamental formula in Refs. [5, 6 does not 
hold as an exact identity and the universal constant reported in Ref. [BJ should be corrected, the logarithmic 
term obtained in Ref. [5J could still stand valid. I will briefly discuss this point in Sec. [BJ) 

For a free boson CFT, new fields defined as linear combinations of the original field are apparently 
independent of each other. However, they are not completely independent since the compactification of 
the new fields intertwine different components. This complication is often ignored in literature and still 
correct results are obtained in some cases. However, its negligence can lead to erroneous results; the present 
problem of entanglement entropy is indeed such an exampleo The intertwining of new free boson fields in 
context of boundary CFT was discussed by Wong and Affleck [5] for a quantum impurity problem. There, 
the compactification of each component of new fields is written explicitly in terms of gluing conditions. 
Although taking all the gluing conditions into account should lead to a correct result, it becomes increasingly 
cumbersome for larger number of components. 

Instead, the compactification of multicomponent bosons can be formulated geometrically using compact- 
ification lattices in multidimensional space. This is useful in construction of some of the possible boundary 



1 The problem in the treatment of compactification was pointed out earlier by V. Pasquier (private communications), as a 
possible source of error in Ref. [6] . 



conditions systematically, and has been applied to string theory and to condensed matter physics. This ap- 
proach is particularly suitable for the present problem, as the replica trick requires calculation for arbitrary 
number of components. 

I demonstrate that the entanglement entropy in the critical groundstate wavefunction of quantum Lifshitz 
univesality class can be indeed calculated fully taking the above subtlety into account, with the geometrical 
construction of the boundary conditions of multicomponent free boson CFT. As a result, I find that the 
universal constant term in the entanglement entropy is in exact agreement with that obtained in Ref. [7] 
by different approaches. This agreement confirms the universality of the constant term, as put forward in 
Ref. |S]. 

This paper is organized as follows. In Sec. [21 1 quickly review the arguments in Refs. 0(3]. I then point out 
subtle problems in the "changing the basis" trick used in these papers. In order to demonstrate the problem 
in the "changing the basis" trick, in Sec. [3l I discuss a simple problem of a single component free boson field 
theory. By a "folding" trick, the problem can be regarded as two-component free boson with boundaries. 
However, naive application of the "changing the basis" technique fails to reproduce the partition function. In 
Sec. 21 I calculate the entanglement entropy for the cylinder geometry, using a geometric formulation of the 
boson compactification. The constant term in the entanglement entropy is given in terms of the universal 
boundary entropy in boundary CFT. The calculation is extended to the torus geometry in Sec.[SJ where it is 
pointed out that the entanglement entropy in the torus geometry is twice of that in the cylinder geometry for 
conformal critical points. Sec. [B] is devoted to conclusion and discussions, including whether there would be 
any correction the logarithmic term predicted in Ref. [5] or not. The recent e-print by Hsu and Fradkin [3] 
will be also discussed in this section. 

Relevant materials in boundary CFT of multicomponent free boson field theory are summarized in Ap- 
pendix [X] In Appendix [B] the solution to the simple example introduced in Sec. [3] is given by explicitly 
solving the gluing conditions. 

2 Setup of the problem 

In a class of two-dimensional critical systems, the groundstate wavefunction is related to a two-dimensional 
CFT. Namely, the ground state is given as 

\^o)^Jv<Pe~ s ^^ 2 \{<P}), (2) 

for the action S[4>] of a CFT. 

I consider a groundstate given as eq. ©. I divide the system into two regions A and B. The entanglement 
entropy between the two regions is defined by the von Neumann entropy of the subsystem A as 

Se = -Tr (p A log p A ) , (3) 

where 

p A = Tr B |* ) (*o| (4) 

I will follow the argument by Fradkin and Moore to derive Se for such a system. [5] Namely, the 
entanglement entropy is rewritten as 

Q dTlpAn I M 

In the replica trick, we compute TvpA n for an integer n and then make an analytic continuation to arbitrary 
integer n. For an integer n, we can introduce n copies of the CFT with the fields </>i, </>2, . . . </>«• 



Trp A n = f^ (6) 



where Z-p is the partition function of the n-component field theory with the condition 

4>l = <j) 2 = ■ ■ ■ = 4>n 



(7) 



at the boundary between A and B. Zp is the partition function of the same n-component field theory but 
without any restriction at the boundary between A and B. Both Zp and Zp are functions of n, although 
the dependence is omitted for brevity of the expressions. 

Fields with different replica indices are independent, except possibly at the boundary T between A and 
B. Since no coupling is introduced at the boundary in Zp, we find 

Z F = {z F ) n , (8) 

where zp is the partition function of the single component free boson field theory without any restriction at 
the boundary. I find no problem in the argument, up to this point. 

They proceed further by changing the basis, taking the linear combinations of the original fields 4>j, as 
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in 



<f X = -L (0x - fa) (10) 

tp n -l = 

It was argued that, the condition does not affect the "center of mass" field <po and it thus remains free 
at r. On the other hand, all the other linear combinations <pj with j > 0, which correspond to differences 
among <fik, obey fixed (Dirichlet) boundary condition at T. The n — 1 "difference" fields was then regarded 
as independent. As a consequence of this argument, it was proposed in Ref. [SJ that 

Z v = (z D ) n ~ z F , (11) 

where zp, is the partition function for the single component field, with Dirichlet boundary condition at the 
boundary. 

Combining eqs. (JS]),©, and (TTTTi leads to 

TrPA " = ( %) ■ (12) 

This implies 

S* = -log2£ = -log^>, (13) 

Zp Zp 

where z^(z^) is the partition function for the single component free boson field theory restricted to the 
region A (or B), with Dirichlet boundary condition at the boundary between A and B. Eqs. (fTTj) and (fT3|) 
are also the basis of the calculations in Ref. [6] . 

However, in a general CFT, the field <fi is subject to interactions. For example, the two-dimensional 
critical Ising model correspond to <f> field theory with a certain fine-tuning. The field theory with the 
interaction is generally not invariant under orthogonal transformations of the fields. For example, let us 
consider the simplest case n = 2. The Lagrangian of two identical Ising field theory would read 

£ (lsing)2 = E (\(Wi? + 5"»V + W) > ( 14 ) 

3=1,2 V - / 

where m 2 and A are fine-tuned to make the system critical. Now I introduce the new fields y>o,i as in eqs. © 
(for n = 2) and (fTU|) . Then, in terms of the new fields, 



£ (ising) 2 = Yl ( o^^ 2 + o m2 ^' 2 ) + I 2 ^° 4 + 6(f ° 2ipi2 + Vl ^ ■ 

.7=1,2 V ' 



(15) 



The new fields ^0,1 are subject to different interaction from the original one. Moreover, different components 
ipo and ipi are now coupled through the bulk interaction ipo 2 ipi 2 and not independent of each other. Thus 
cq. pip would not hold in general field theories with interactions. 

For a free field theory, on the other hand, the change of the basis appears more legitimate, because the 
theory is also free in terms of the new basis fields. Let us now consider the case of the free boson field theory. 
As I have discussed, this will be relevant for entanglement entropy in the quantum Lifshitz universality class. 
I define the Lagrangian density following the convention in Ref. [10] as 

£=£(^) 2 - (16) 

The field cf> is subject to compactification, namely the identification 

<f> ~ <f> + 27ri?, (17) 

where R is the compactification radius. 

I note that, there is unfortunately a large variety of conventions for the free boson field theory, based on 
different normalizations. In fact, by a renormalization of the field <f>, we can fix the value of either coupling 
constant g or the compactification radius R. However, we cannot fix both g and R by the renormalization. 
This leaves one free parameter, which governs the critical behavior. We can choose either to fix g and 
consider R as a free parameter, or to fix R and regard g as a free parameter. Both conventions, with various 
choices of the fixed value, appear in literature. In this paper, I keep both g and R as parameters. This is 
a redundant parametrization, but makes it easier to compare with literature by setting either g or R to the 
convention of each paper. For example, the convention in Refs. [SI [7] can be recovered by setting g — 1/2. 

Employing the replica trick, I consider n component free boson theory, in which each component is 
independent in the bulk. For the free field case, the Lagrangian density can be also written as 

n-1 
£ =£l>/^) 2 - ( 18 ) 

That is, the new fields are governed by the same Lagrangian density as the original fields, and there is 
no interaction which couples different components. Thus the arguments in Refs. [SJ [5] appears to valid. 
However, even for the free boson field theory, eq. (jlip docs not quite hold because of subtlety in boson 
compactification. 

3 A simple example 

In order to illustrate the issue, let me discuss a simple example. I consider a single component free boson 
theory (fT6]) on a rectangle of the size 2L x j3, with the periodic boundary condition on both directions. In 
other words, the system is defined on a torus. The partition function is given by 

-^simple — /%,_% / <i H 

■q(q) T]{q) ^ 



' \\^ +n2 ° R2 \ (19) 



v(v) 



where 



-77/3/L 



(20) 

For a general torus, q is given as q — e 27rir where r is the modulus of the torus which is a complex number, 
and its complex conjugate q is distinguished from q. However, here I only consider the rectangular case 
where r is a pure imaginary, and hence q — q. I have used this fact in the second line of eq. (|19l) . 




Figure 1: A single component free boson field theory is defined on the torus of size 2L x f3. The system is 
folded to a two-component free boson field theory defined on the cylinder with circumference /3 and length 
L. 



Now I "fold" the system, as shown in Fig. Q] After the folding, the system can be regarded as two- 
component boson field 4>i,4>2 defined on a rectangle of size Lx/f. While the periodic boundary condition 
is still applied on the /3 direction, there are two boundaries which were the folding lines. Thus the system 
after the folding is topologically a cylinder. At the two boundaries, I impose the condition 



(21) 



Namely, if we consider the new fields 



*' = ^TT- (23) 

<&i obeys the Dirichlet boundary condition $1 = 0, while $o remains free (obeys the Neumann boundary 
condition) at the two boundaries. 

If we apply the same argument as in Ref . [6] , we would obtain 

Zsimpic^) = z DD (R,q)z NN (R,q), (24) 

where 

z DD (R,q) = 4r £«•"'* = ^T^E^ 49 ^ (25) 

is the Dirichlet-Dirichlct amplitude for the single component boson with the compactification radius R, and 

is the Neumann-Neumann amplitude for the same theory. Here, 77(g) the Dedekind eta function defined in 
eq. (|98|) . However, from eqs. (|25|) and (l26t we can immediately see that eq. (l24l actually does not hold. 
One of the problems is that, in the "closed string channel", the right-hand side of eq. (|24|) reads 

1 / 1 ^ 2 



znn(R,q)z NN (R, q ) = ^[^) E«" '^ )+ ™ ^ " ' ^ 



li.rn 



This implies that each boundary has the "groundstate degeneracy" (exponential of the boundary entropy) 
l/v2, for any value of R. This is in contradiction to the fact that the boundary in the present example 
is just a result of an artificial "folding" along a line in the bulk, and thus should not have any boundary 
entropy. In fact, the modular invariance of the partition function on the original torus of 2L x ft implies 

^-(^Ef^^'- (28) 

with the coefficient unity. Interpreted as the amplitude with two boundaries after the folding, this means 
that the boundary entropy is indeed zero. Thus, there must be something wrong in the assumptions led to 
eq. (|24p (and to eq. (|lll) in Ref. [B]) even for the free field. 

The problem was that the new fields were implicitly assumed to obey the same compactification as the 
original fields: 

$j - $, + 2irR. (29) 

In fact, while the original fields arc compactified independently, a complication is introduced in the com- 
pactification by a change of the basis. [5] This can been seen in Fig. [21 The compactification in terms of the 




Figure 2: The compactification lattice for two independent boson fields 4>i,2, each of which has compactifi- 
cation radius R. Points on the compactification lattice (shown as red and blue circles), which is the square 
lattice with lattice constant 2irR in the (<f>i, 02)-plane, are identified. The compactification is not imposed 
independently on the linear combinations $o = (4>i + 4>2)/V^ and $1 = (cf>\ — (j>2)/^/2. The compactification 
lattice can be divided into two sublattices (red and blue). Considering only the red sublattice would be 
equivalent to independent compactification of $o and $1 with radius \/2R. However, the blue sublattice 
should also be included in the compactification lattice. 



new fields $0,1 reads 



where 



$o~$o + 27m ^=, (30) 

$i~$i + 27rni^, (31) 

v2 



no = n% mod 2. (32) 



Here no = n i = mod 2 and no = n\ = 1 mod 2 correspond respectively to the red and blue sublattice in 

Fig.m 

We can also define the new fields 80,1 similarly for the dual fields 9\^- Similarly to the case of $0,1; their 
compactification is given as 

e ~e + 27rmo-=— , (33) 

91-91 + 27,-™!—!—, (34) 

V2gR 

where 

ttt-o — m i mod 2. (35) 

These considerations imply that the compactification is not independent in terms of the new fields and 
is subject to "gluing conditions" (|32t and (|35|) among different fields. This aspect was ignored in Ref. [6j. In 
the present simple example, we can see that ignoring the gluing conditions leads to a wrong equality (|24|) . 

This simple example demonstrates the importance of the gluing conditions - namely, that the linear 
combinations of the compactificd fields are not completely independent. In fact, taking the gluing conditions 
into account, the correct partition function can be reproduced from the boundary CFT. I describe this 
calculation in Appendix [Bl 

In passing, I note that in Ref. [8] additional gluing conditions between the dual winding numbers (rij and 
rrij) were discussed. These originate from Fermi statistics of electrons in the microscopic model. (See also 
Ref. [ID].) Those extra gluing conditions do not apply to the present case. 

4 Cylinder geometry 

Now let us move on to the problem of the entanglement entropy. I start from the "cylinder" geometry 
introduced in Ref. |6 , as shown in Fig. |3] Namely, I consider a cylinder of circumference j3 and length 
La + Lb- It is divided into two regions A and B, with length La and Lb respectively. For simplicity, I 
consider the case La = Lb = L. At the two ends of the cylinder, I impose the Dirichlet boundary condition 
4> = 0, as in Ref. [6j. 

Following Refs. |6j, I employ replica trick. At the boundary T between A and B, the condition ([7]) 
is imposed. In order to apply boundary CFT to the present situation, I invoke the folding technique as 
introduced in Ref. [H] and also discussed in Sec. El Then the system may be regarded as a 2n-component 
free boson field theory on a cylinder of length L and circumference (3. I simply label the fields after the 
folding as <fij, where j = 1,2,..., 2n. Here <j)j for j > n represents the "folding double" of (j>j- n - I will denote 
the doubled number of components 2n as J\f. Each component obeys the compactification as in eq. (|17|) . 
Thus, the compactification lattice A is the A/"-dimensional hypercubic lattice with the lattice constant R. 

At one end of the cylinder, which corresponds to the two ends of the original cylinder before the folding, 
the Dirichlet boundary condition 

(j>j = (36) 

is imposed for j = 1, 2, . . . , N . 




Dirichlet b.c. 



r 



A / 




B 




-►-«- 



Dirichlet 
b.c. 




folding 





r 



Figure 3: The upper panel shows a cylinder of circumference [3 and length 2L. Dirichlet boundary condition 
<fi = is imposed at the both ends. The system is divided into two regions A and B with length L each, and 
I discuss the entanglement entropy between the regions A and B. In replica trick calculation, n-component 
free boson field theory is defined on the cylinder. It is folded onto a cylinder of length L, with the open 
boundary condition on the one end. The other end corresponds to the boundary T between the regions A 
and B. After the folding, M = 2n component free boson is defined on the cylinder. 
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In calculating Z-p, the condition (J?} is imposed at the other end, which corresponds to the boundary 
between the two regions. In terms of the A/"-componcnt field, it reads 

<fa = 02 = • • • = 4>M-i = 4>M- (37) 

While j\f = 2n is an even integer for the present application, the following construction of the boundary 
state is valid for any (positive) integer j\f. 

Following Refs. 0E], we ma y define the new basis by 



i 

$n = 



3 = 1 

*isi (0i ~ 02) (39) 



$A/--1 = • • • ■ 

We can also define their dual, Oo, . . . , 0jv"-i by linear combinations of 9j. Among the new fields $j, only $o 
obeys free (Neumann) boundary condition, while all the others $i, $2, • • • $jV-i is subject to the Dirichlet 
boundary condition <frj = 0. Such a "mixed" (Dirichlet/Neumann) boundary condition has been discussed 
in string theory [12] and in condensed matter [10] applications. 

As we have seen in Sec. [3] (which corresponds to Af = 2), the compactification is not independent in terms 
of new fields $j. It follows that the "mixed" boundary condition cannot be given just by a simple product 
of Dirichlet/Neumann boundary states. For a correct evaluation of the entanglement entropy by the replica 
trick, we need to take gluing conditions among Af = 2n component bosons correctly into account. However, 
it is cumbersome to keep track of the gluing conditions explicitly as in Appendix [B] for larger number of 
fields. Fortunately, the compactification of multi-component boson field can be handled systematically with 
a geometric formulation of the "compactification lattice" . 

I will denote the boundary state corresponding to this "replica" boundary condition as \V). The boundary 
condition implies that, the orthogonal matrix 72 in eq. (|82|) is reflection about the plane normal to the j\T- 
dimensional vector 

d=(l,l,l,...,l) T . (40) 

Explicitly, the matrix is given as 

72 = 1 - 2d(F. (41) 



In order to construct the boundary state, we need to identify the winding numbers which satisfy eq. 
I define 'B.dv as the intersection of A* / \fg and the M — 1 dimensional hyperplane d ■ K = 0. It is a Bravais 
lattice on the hyperplane d ■ K — 0. The general solution (R, K) of eq. (|86p. for the current choice of 72 in 
eq. (|4Tj) . is given by 

R = n Rd, (42) 

for any integer no, and any K € H^p. 

The corresponding boundary state may be written as 

\V{B Q M))=9v E e-^-^^ I (£,#))>, (43) 

R— noRd,K£'ED-p 

where the summation is taken over the solution of eq. (|86p discussed above. 9q and cj>o corresponds to the 
boundary values of 9 and </>. Since R \\ d, components of #0 which are orthogonal to d is irrelevant in this 
expression. Thus we can assume without losing generality that 9q is parallel to d. Similarly, components 
of 4>o parallel to d is irrelevant and thus <fio can be assumed to be orthogonal to d. These correspond to 

11 



the fact that the present boundary condition is Neumann on $o (i*e. Dirichlet on Oo) and Dirichlet on 
$1, $2, • • • , $jV-i- $o and </>o represent the boundary value of Oo and $i t ... ,.aa-i respectively. 

Actually, the "replica" boundary condition implies all the "difference" fields 'E'i,... ,7V— l vanish at the 
boundary, and thus all components of 4>q is zero for the present problem. However it is useful to remember 
that this boundary state is a special point in the continuous family of boundary states labelled by 9q || d and 
4>q _L d. This demonstrates that \P) is the boundary state indeed corresponding to the Dirichlet boundary 
condition for Go and $i,2,...,jV-i- 

The coefficient g-p is determined by Cardy's consistency condition. The numerator of eq. ^ in this case 

is given by the amplitude of the A/"-component free boson field theory with the boundary states \D) and \P) 

at the two ends. Only the winding number sectors common to \D) and \V) contribute to the amplitude. 

Thus 1 find 

/ i \ ^~ x °° 1 

m — l iieSn-p 

In order to satisfy Cardy's consistency condition upon the modular transformation, the amplitude must be 
written as 

Zdv{q) = z DN (q)f(J\f - 1; E DV ;q), (45) 

where J- is defined in eq. (J96I) . By comparison to eq. (|96[) . I obtain 

g D g v = ^-^-^voiZn-p) (46) 



In the present case, 
which implies 



g D = {2g)- M l 4 R- M ' 2 (47) 

9v = (^^R^VoPdv) (48) 

Now I determine Vo{Zdv)> by calculating w (A=) in two ways. First, since A* is the A/"-dimensional 



hypercubic lattice with the lattice constant l/R, we find 

A* / 1 \ M 

On the other hand, vq(A=) can be also written in terms of voi^Dv), as follows. The Af — 1 dimensional 

lattice "E.DV is a intersection of A* and the hyperplane orthogonal to d. This is illustrated in Fig. [4] for the 
simplest case of N = 2. The hypercubic lattice A* / \Jg can be decomposed into parallel displacements of 
"B.DV- Because A* is the simple cubic lattice with lattice constant 1/R, for any vector K e A* /-Jg, 

-dm , 



with an integer m. On the other hand, for any integer to, there is always a vector K £ A* / \fg which satisfies 
eq. ([5U| . Eq. (|5U| represents the distance between the hyperplane and the origin. When eq. (fST)]) holds, K 
belongs to the m-th parallel displacement of Zdv ■ Thus the distance between the neighboring hyperplanes 
hosting a copy of S^-p is 

7m- <51) 

and the volume of the unit cell of the original lattice A* /.Jg is given as 

v (^=) = —L-v (E DV ). (52) 

yfg y/gJVR 
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^DV 



Figure 4: The construction of the lattice Eov is shown for J\f = 2. The lattice A* /<Jg is A/"-dimensional 
hypercubic lattice (square lattice for J\f = 2) with lattice constant 1/^/gR, shown by red and blue circles. 
'E.dv is given by the section of A* /^/g by the Af — 1 dimensional hyperplane (line for N = 2) orthogonal 
to d = (1,1,..., 1) T , shown by a yellow line. A* /^/g can be decomposed into copies of Eov displaced in 
parallel. The distance between the neighboring copies of Edv is given by l/(\/gAI B). 
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By comparison of cqs. (l4"9")l and ([52")l . I find 

v (~ DV ) = y/2g-W-WRTW-»yJ¥ (53) 

Therefore eq. (|48fl implies 

gp-^Kl-U- 1 ^ (54) 

= (v^iJ)"^ 1 ^ (55) 

The denominator of eq. (JSJ corresponds to the amplitude Zuf of the A/"-component free boson field theory 
with the boundary states \D) and \F) at the two ends. \F) is the boundary state corresponding to 



"j+rn 



(56) 



for j = 1, 2, . . . , n. Namely, each field is identified with its "folding double" at the boundary. This corresponds 
to a boundary created by folding n-component free boson in the bulk without any defect. 

As I have discussed in Sec. |31 the boundary entropy for this artificially created "boundary" should be 
zero for any n. Thus it must follow that 

9F = 1, (57) 

for any n. In fact, eq. (|57|) for the boundary state \F) can be also shown by an explicit calculation similar 
to that of g-p. I note that, both \F) and \V) reduces to \P) for n = 1 (W = 2). 
In the long cylinder limit L 3> /3, q — > and I obtain 

z v = z DV (q) - q- ,l/12 g D gv, (58) 

z F = z DF (q) ~ q- n/12 g D gF. (59) 

Thus, using eq. ([57]). I find 

Tr PA "~ 5P , (60) 

in the long cylinder limit. The von Neumann entanglement entropy from the replica trick is thus given as 

From eq. (1551) . for the quantum Lifshitz universality class, 

S E = log (y/2JJR) - \. (62) 

To match the convention in Refs. [6l [7], I take g = 1/2 and thus Se — log J? — 1/2. This is different from 
the result reported in Ref. [6], by the second term —1/2. On the other hand, it indeed agrees exactly with 
that in Ref. 7 derived by different approaches. This shows that the basic ideas put forward in Refs. [51 [3] 
are correct, although care must be taken in changing the basis. Several results reported in Ref. [3] are not 
valid owing to negligence of the subtlety. The validity of the logarithmic term predicted in Ref. [5] will be 
discussed in Sec. [SI 

5 Torus geometry 

In Ref. [5], the entanglement entropy in the torus geometry was also discussed. Namely, a torus of total 
length La + Lb is divided into two regions A and B of cylindrical shape. The boundary T between A and 
B consists of two disjoint circles. 
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L 




folding 




Figure 5: The torus (upper panel) of size j3 x 2L is divided into two regions A and B with length L each, and 
I discuss the entanglement entropy between the regions A and B. Folding reduces the torus to the cylinder 
of length L, where the boundaries at the two ends correspond to the boundary T between regions A and B. 
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Again, for simplicity I consider the case La = Lb = L and apply the folding trick to reduce the problem 
to JV — 2n-component free boson field theory on a cylinder of length L, as in Fig. [5] By construction, in 
the present case, the "replica" boundary condition V should be imposed on the both boundaries, in order to 
calculate the numerator of eq. ©. This is contrasted to the case of cylinder geometry discussed in Sec. |H 
where the replica boundary condition is imposed at one end and the Dirichlet boundary condition at the 
other. 



Thus I find 
In the long cylinder limit, 

This implies that 
which leads to 



TvpA = zm$- (63) 



Trp A »~(?Z) =g v \ (64) 

g(torus) _ 2g(cylindcr) /ggx 



4 tOTUs) = 2 (lo gv /2^-ij. (66) 

This result is, again, different from that in Ref. [5J. In particular eq. (p5|) does not hold in Ref. [BJ. On the 
other hand, eq. (l6"5l) is a general consequence of the boundary CFT, independent of the value of Se- In fact, 
although it was not explicitly discussed in Ref. [7] , eq. (|65l) is also a general consequence of the mapping to 
classical statistical problem used in Ref. [7] . Thus the violation of the relation (|65|) is a clear signature of 
the problem in the calculation in Ref. [BJ. 

I emphasize that eq. (1651) does not hold for the entanglement entropy in general systems. It only applies 
to the critical wavefunction described as in eq. ([%]). for which the entanglement entropy can be related to 
classical (Shannon) entropy. 7 In fact, eq. (|65[) is violated in topological Z2 spin liquid phase. 

6 Conclusion and Discussions 

I have discussed the entanglement entropy in two-dimensional conformal critical points, in particular those 
described by free boson CFT (quantum Lifshitz universality class). 

Calculations in Refs. [5j[6j were based on the fundamental formula ([TT]) (or equivalently eq. ((T3j)) . How- 
ever, there are two problems in this formula. First, in a general interacting theory, linear combinations of 
the original fields are not independent of each other as it was implicitly assumed in Refs. 5,6). On the other 
hand, CFTs other than the free boson field theory also admit free field representations. It might be used 
to extend the approach to general CFTs. However, it would be a nontrivial problem which requires further 
careful investigations. 

Second, for the free boson field theory (which describes the quantum Lifshitz universality class), the bulk 
interaction is absent and the linear combinations appear independent. However, even in this case, they are 
not completely independent because of the intertwined compactification. I have demonstrated the importance 
of the nontrivial compactification using the simple example of single-component free boson field on a torus, 
which can be regarded as two-component free boson field with boundaries. "Mixed" Dirichlet/Neumann 
boundary conditions for general number of components, which appears in the replica trick calculation of the 
entanglement entropy, can be handled with a geometric formulation based on the compactification lattice. 
The constant part in the entanglement entropy corresponds to the universal boundary entropy in the CFT. 
This supports the universality of the constant term as proposed in Ref. [6j. In fact, the entanglement 
entropy obtained for the quantum Lifshitz universality class with the boundary CFT agrees exactly with 
that obtained with different methods in Ref. [7] . An incorrect value was reported in Ref. [B] because of the 
compactification was not properly taken into account. 

The predictions made in Ref. [5] on the logarithmic term in eq. (U} should also be re-examined, since the 
fundamental formula dill) , on which their derivation is based, does not hold as an exact identity. Nevertheless, 
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their prediction on the logarithmic term could still stand valid for the following reason, in particular for the 
free boson (quantum Lifshitz) case. For the free boson field theory, in Ref. [5 , they found the logarithmic 
term independent of the compactification radius. This suggests that the logarithmic term could be attributed 
solely to oscillator modes. For oscillator mode contributions, there is no subtlety discussed in the present 
paper and the "changing the basis" trick could be justified. From this perspective, it seems quite possible 
that their prediction on the logarithmic term is correct despite the subtle problem with eq. dill) . It would 
mean that eq. (Ilip is valid in some restricted sense for determining the logarithmic term, although it is 
certainly not valid for determining the universal constant term. For general CFTs other than free boson, 
there is more problem in the "changing the basis" trick due to the bulk interaction. It might be still possible 
that the logarithmic term depends only on the central charge and the prediction in Ref. [5] is also valid for 
general CFTs, but it seems less convincing than in the free boson case, at this point. 

In any case, in the present paper, I do not have any concrete result concerning the logarithmic term, and 
cannot draw a definitive conclusion about the prediction in Ref. [5] on the logarithmic term. I hope that the 
present paper will stimulate further progress in understanding of the logarithmic term. 

It should be noted that the present paper, as well as the original work [B], is entirely based on (a simple 
implementation of) the replica trick. Its validity is by no means obvious. For the present case of free 
boson CFT, the agreement with different methods [7J implies that it is indeed valid. However, recently, its 
breakdown is suggested when the CFT corresponds to critical Ising model. [13] A general understanding of 
the issue is an important open problem. 

Finally, when this paper was close to completion, a paper by Hsu and Fradkin has appeared [5]. There, 
they did not rely on eq. (ITBl . which is not valid as I have discussed. Instead, they attempted to construct 
the boundary state in the spirit similar to the present paper. However, although their construction (eq. (14), 
or eq. (18) with eq. (19) of Ref. [9]) should give a consistent boundary state, it does not correspond to the 
required boundary condition (|37j) for the problem. As an indication, the condition (|37|l is invariant under 
any permutation of boson fields (f>j while their construction is not. 

Furthermore, apparently there are several confusions in Ref. [9]. For example, in the original paper [6], as 
well as in Ref. [7] and in the present paper, the universal constant in the entanglement entropy was derived 
in the "long cylinder limit" L ^> (3 (in the notation of the present paper). However, in Ref. [BJ, the definition 
of the lengths is somehow exchanged and the opposite limit is taken. As discussed in the present paper, in 
the long cylinder limit as introduced originally in Ref. [BJ , the universal constant part of the entanglement 
entropy should correspond to the universal boundary entropy (exponential of the groundstate degeneracy) 
in the boundary CFT. 
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A Boundary CFT of multicomponent boson field 

In this Appendix, I summarize the relevant formulae in the boundary CFT of multicomponent compactified 
free boson field theory. 

The boundary CFT was largely developed by Cardy Q31 [15] . The development of boundary CFT of free 
boson field theory was also started in the context of string theory [HI [17] , and continued for example in 
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Refs. [TSJ dH HOI HI] ■ For a review in the string theory context, see Refs. [2"T] . 

On the other hand, the boundary CFT is also relevant for some problems in condensed matter physics, 
such as quantum impurity problems and junction of one-dimensional quantum systems. The boundary 
CFT of free boson field theory is applied to impurities in quantum spin chains by Eggert and Affleck 22 . 
The boundary CFT of multicomponent free boson is applied to impurities in quantum wires by Wong and 
Affleck [5] , and further discussed in related problems [23, |2H [TU] 

The multicomponent free boson field theory is defined by the Lagrangian density 



where <f> is a A/"-dimensional vector. 

I introduce the multidimensional generalization of the boson compactification as 

6~$+2ttR, 



(67) 



(68) 



where R e A for a Bravais lattice A, which is called as compactification lattice. For Af independent copies of 
the single-component boson compactified as in eq. (|17|) , the compactification lattice A is simply a hypercubic 
lattice with lattice constant R. However, it is useful to formulate allowing more general compactification 
lattice. 

I note that, as in the single-component case, the coupling constant g can be set to any value by renor- 
malizing the field 0, which also renormalize the compactification lattice. Thus g is a redundant parameter 
once I consider general compactification lattice. Nevertheless, I keep the coupling constant g for purpose of 
comparison as I have discussed for the single-component case. On the other hand, for general compactifica- 
tions, it is impossible to fix the compactification lattice by renormalizing <f> within the Lagrangian density of 
the form (|67|). 

If we define the theory (|67|) on a finite length fi with the periodic boundary condition, the canonically 
quantized operator (f> is given by the mode expansion 



0(t,x)=^> + — 



2n 
J 



Rx + Pt 



1 



1 



2<? ^i 



a„exp 



2tt 
T 



a^exp 



2tt 
J 



(69) 



where x± = x±t. R G A represents the winding number of <f> when x goes around the system. The canonical 
commutation relation of <f> with the conjugate momentum field implies 



bf\P k ]=iS jk . 



(70) 



Since the constant part 4> is a l so subject to the compactification as in eq. (|B5|) . the eigenvalues of the 
conjugate "momentum" operator P is quantized as 



9 



(71) 



where K belongs to the dual lattice A*, which is defined as a set of all Af's which satisfies 

K ■ R — integer, 



(72) 



for any R £ A. The momentum operator P in eq. (|69j) is often rewritten using eq. (|7I[) . However, it should 
be kept in mind that K represents the eigenvalue of P as in eq. (|7"Tj) . 

The boson field <fi can be decomposed into left-moving and right-moving components as 



'(x+) + <f, R (x-). 



(73) 
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We can introduce the dual boson field as 
From eq. (I69|) , mode expansion of 9 is given as 



e{t,x) = m + 2 j f 



Kx + gRt 



'IE4 « 



exp 



2tt 
-znx+ — 

P 



a R exp 



2tt 
J, 



(74) 



(75) 



In this expression, the roles of K £ A* and gR € <?A are interchanged compared to the mode expansion of < 
in eq. (|M)l . This implies that the dual field obeys compactihcation 



. 2tt - 
+ —K, 
9 



(76) 



where K G A*, gi? is now interpreted as the eigenvalue of the dual "momentum" operator, which obeys 
canonical commutation relation similar to eq. (|70p with the operator 9^°' . 
The Lagrangian density (|57| can be also written in terms of 9 as 



1 
47rg' 



/: = — -{djf. 



(77) 



Although 4> and are mutually non-local, construction of a complete set of physical operators requires both 
</> and 9. 

When we take boundary of a 1 + 1 dimensional field theory orthogonal to the (imaginary) time axis, the 
boundary condition can be represented by the initial (or final) state. This is called as boundary state and 
is useful in systematic study of boundary conditions. In fact, nontrivial boundary conditions often cannot 
be defined precisely without introducing the boundary states. When the boundary condition is conformally 
invariant, the corresponding boundary state B) satisfies 



L, n 



L_ 



l-B) = 0, 



(78) 



for any integer m, where L m are generators of the Virasoro algebra. For the multicomponent boson field, 
the Virasoro generators are given as 



iJyn. — 



l m-l^l 



-Y, 

i 



(79) 
(80) 



where 



ok = 



—iy/nd^ (n > 0) 




j ~iVna% 


(n > 0) 


75(^+71*0 ^ = °) 


, <** = 


75 (-v^+7^) 


(n = 0) 


i^{a L _J (n < 0) 




{ Wni&nV 


(n<0) 



(81) 



Here each component of a^ ,R is a boson annihilation operator, RgA, and K e A*. a^ R and their Hermitian 
conjugates (boson creation operators) represent the oscillator modes, corresponding to the quantized normal 
mode oscillations. These oscillator modes do not depend on the compactification. On the other hand, a L ' R 
represent zero modes which are affected by the compactification. 

The general solution of the conformal invariant boundary state (|78[) for the multicomponent free boson 
is not known. However, a sufficient condition for eq. (|78[) can be given as 



(at - nd R 



\B) = 0, 



(82) 
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for arbitrary integer m. Here 72. is an A/" x A/" orthogonal matrix independent of m. Eq. ([52")) for m 7^ 
determines the boundary state to be of the form 

expUf>£)t72(a^Avac>, (83) 

where |vac) is an oscillator vacuum. This is a free boson version of Ishibashi state |25j . which is conformal 
invariant. 

In fact, there is an infinite number of oscillator vacua characterized by the zero mode quantum numbers. 
Following Ref. [10 , I label these vacua as 

\(R,K)). (84) 

R and K can be interpreted as the "winding numbers" of (p and of its dual 6, respectively, along the boundary. 
The Ishibashi state obtained as eq. (1531 from eq. (|84"]) is denoted as 



\(R,K))). (85) 

Eq. (|82[) for n — puts a restriction on the vacua which can appear in the boundary state. That is, the 
quantum numbers must satisfy 

^R+^k\=Tz(-^R+^-k\. (86) 

For given orthogonal matrix 72, generally there is an infinite number of solutions (i?, K) which satisfy this 
requirement. Any linear combination of the Ishibashi states built from these vacua satisfies the conformal 
invariance (175)) . 

However, a physical boundary state must also satisfy Cardy's consistency condition, stated as follows. 
For a pair of given boundary conditions A and B, we can define the amplitude (partition function) 

Z AB {q) = (A\e- Lfl -\B) (87) 

where q — ... By modular transformation, we can express this amplitude as a function of q = e, 

ZAB(q)=Y, N ABxf r (q) (88) 

h 

where Xh^il) i s a character of the Virasoro algebra. Since Nhg can be interpreted as the number of primary 
fields with conformal weight h, it has to be a non-negative integer. This is Cardy's consistency condition. 
Usually it is also required that N^ . — 1, where h = corresponds to the identity operator. If one takes 
just a single Ishibashi state, Cardy's condition cannot be satisfied. Generally, we take linear combination 
of Ishibashi states for all the zero modes allowed by eq. (1561 . It is also possible to construct consistent 
boundary states using only subset of zero modes allowed by eq. (|5B"|) . However, those boundary states are 
more unstable and the most stable boundary states for a given 72 turns out to be linear combinations of all 
the allowed zero mode vacua [TU] , 

For example, taking 72 = 1 gives Dirichlet boundary state. The solution of eq. (156")) for 72 = 1 is given by 
R = 0. Thus the Dirichlet boundary state is given as 

\D($ ))= 9D J2 e- t$ °- R \(0,K))), (89) 

where the summation in K is taken over the entire dual lattice A*, and 4>o is a constant A/"-dimensional 
vector. (There is an unfortunate conflict of notation; go here is the coefficient of the boundary state and is 
a completely different quantity from the coupling constant g defined in eq. Ijl6[) .) 
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Physically, the constant vector </>o corresponds to the boundary value of the field <j). This can be seen as 
follows. Let us consider the operator 

where K$ is a constant vector belonging to the dual lattice A* . This is a unique valued operator under the 
compactification[BSJ We apply this operator to the boundary state ([55)) at (imaginary) time 0. With respect 
to oscillator modes, each Ishibashi state \(R,K))) constructed as in eq. (|83[) is a kind of coherent state. We 
observe that, for general H, 

a L n \{R,K)))=-n{a^\(R,K))). (91) 

For the Dirichlet boundary condition, 1Z = 1 and thus applying the oscillator part of 4> to the Ishibashi state 
yields zero. The winding number R is zero in all the Ishibashi states in the boundary state (|59")l . and at 
t = 0, contribution of the "momentum" eigenvalue K to vanishes. Thus we find 

e i${ - x ' 0) - Ro \D{$ Q )) =e l$m - Ro \D($ Q )). (92) 

Now the commutation relation (1701) implies 

e^ {0) -^\(R : K)) = \{R,K + K )). (93) 

Combining this with eq. (|92[) , we find the eigenequation 

e<#*.°>-*o (£)(£,)) = e l ^ Ro \D($ )). (94) 

This proves that (/>o can indeed be interpreted as the boundary value of the field cf>. The Dirichlet boundary 
state is actually a continuous family of boundary states parametrized by 4>o- By symmetry, physical properties 
such as boundary entropy and scaling dimensions of the boundary operators are independent of (f>Q. Thus 
the boundary value is often set to zero for simplicity. I denote \D((fio — 0)) simply by \D). 
The Dirichlet-Dirichlet amplitude is given as 

Z DD (q)=T(N;^-;q), (95) 

where 

HM-~,q) ^ 2^/ 2 , (S) ( -L Yj2^ /4 (96) 

^ <Af;H ' ;,) ^(^))'£' s "- (97) 

Here S is a Bravais lattice and S* is its dual, i>o(S) is the volume of the unit cell of the lattice S, and 

r)(q)=q 1/2i li(l-q n ), (98) 

is the Dedekind eta function. I note that 

v (E*) = -^-. (99) 

Equality between eq. (|M|) and eq. (|9"Tf holds thanks to the multidimensional generalization of Poisson 
summation formula. 

The coefficient in eq. (|96| is fixed so that the amplitude in the form (|97|) satisfies Cardy's consistency 
condition. 
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By comparison to cq. (|96[) . we find 

g D = (2gy^/\vo(A))- 1 / 2 , (100) 

where i>o(A) is the volume of the unit cell of the compactification lattice A. go can be interpreted as 
universal non-integer groundstate degeneracy; in other words, log go is the boundary entropy for the Dirichlet 
boundary state. 

Likewise, the Neumann boundary condition corresponds to 1Z = —1. The Neumann boundary condition 
for <j) is equivalent to the Dirichlet boundary condition for the dual held 9. The Neumann boundary state is 
given as 

\N(9o))=g N J2 e- i2 ^°-*\M)), (101) 

ReA 

where 9q is the boundary value of the dual field 9. The Neumann-Neumann amplitude reads 

Z NN (q)=f(Af;^gA;q), (102) 

It follows that the groundstate degeneracy is 

g N = {\Y'\v {K)f'\ (103) 

The Dirichlet-Neumann amplitude is also of interest; it must satisfy Cardy's consistency condition as well. 
I note that, the Hamiltonian time evolution does not change the winding numbers from those in the initial 
state. Thus the only oscillator vacuum which contributes to the amplitude is |(0, 0)), and no summation over 
zero modes appears in the Dirichlet-Neumann amplitude. The amplitude is thus given only by oscillator 
mode contributions as 

Z ND (q) = (z ND (q)f , (104) 

where 

= 24)«'" 2 >' < io5 » 

is the Dirichlet-Neumann amplitude for a single component boson. Here $2,4 are Jacobi's theta function 
defined as 

00 

(106) 

(107) 

This amplitude indeed satisfies Cardy's consistency condition, as it should. It is also noted that this amplitude 
does not depend on the compactification. 

B Construction of the boundary state for the simple example 
based on the gluing condition 

Here I construct the boundary state \P) for the boundary condition (|2ip. using explicitly the gluing condi- 
tions ([3"2l and (f3"5j) . Although the boundary state can be given by a more systematic geometric formulation as 
in Sec. the present construction would be instructive to understand the importance of the gluing conditions. 
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V^ q (n+l/2) 




71=— OO 




OO 


^4 (<Z) 5 


E (- 1 ) n « n: 



As I have discussed in Sec. [31 the relevant boundary condition would correspond to the Dirichlet boundary 
condition for $1, namely $1 = 0, and the Neumann boundary condition for $o- The latter is equivalent 
to the Dirichlet boundary condition for the dual field, Go = 0. Following the standard construction of the 
Dirichlet/Neumann boundary state for a free boson field theory, I can write down an ansatz: 

\P)=9p Y, l(no.0,0,mi))), (108) 

no, mi 

where |(?loj n i-> m o? m i)) is the Ishibashi state with the winding numbers along the boundary 

A$ j = 2nnj- F , (109) 

v2 

Ae - 2 ™'db' (110) 



where j = 0,1. The gluing conditions (|32|) and f|35[) imply, because n\ = ttiq = 0, that no and mi must be 
even integers. Thus the partition function (the amplitude with \P) boundary state at the both ends) reads 



Modular transforming, I obtain 



Z PP = ( 9P f(-±-r) Y^ 9Rn2 + ^\ (H2) 



which implies gp — 1 due to the Cardy's consistency condition. Namely, there is no "boundary entropy" 
for this boundary, as it is required by physical grounds. Moreover, the partition function Zpp indeed agrees 
exactly with the original expression (|19j) . 
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